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£j ■ Abstract 
O 

We consider a two-dimensional random resistor network (RRN) in the pres- 

; 

ence of two competing biased percolations consisting of the breaking and re- 

a, 

, covering of elementary resistors. These two processes are driven by the joint 

C , 

effects of an electrical bias and of the heat exchange with a thermal bath 

O 



The electrical bias is set up by applying a constant voltage or, alternatively 
a constant current. Monte Carlo simulations are performed to analyze the 
network evolution in the full range of bias values. Depending on the bias 
strength, electrical failure or steady state are achieved. Here we investigate 
the steady-state of the RRN focusing on the properties of the non-Ohmic 
regime. In constant voltage conditions, a scaling relation is found between 
< R > I < R >o and V/Vq, where < R > is the average network resis- 
tance, < R >o the linear regime resistance and Vo the threshold value for the 
onset of nonlinearity. A similar relation is found in constant current condi- 
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tions. The relative variance of resistance fluctuations also exhibits a strong 

nonlinearity whose properties are investigated. The power spectral density of 

resistance fluctuations presents a Lorentzian spectrum and the amplitude of 

fluctuations shows a significant non-Gaussian behavior in the pre-breakdown 

region. These results compare well with electrical breakdown measurements 

in thin films of composites and of other conducting materials. 
PACS: 77.22.Jp; 64.60.Ak; 07.50.Hp; 64.60.Fr 
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The study of electrical and mechanical stability of disordered systems is attracting a 
considerable interest in the recent literatureEhU because of its implications on both ma- 
terial technoloe and fundamental aspects related to the response of these sys- 
tems to high external stresses@~0'Eil _ i!. Indeed, the application of a finite stress (elec- 
trical or mechanical) to a disordered material generally implies a nonlinear response, 
which ultimately leads to an irreversible breakdown (catastrophic behavior) in the high 
stress limit0~i'0. Such catastrophic phenomena have been successfully studied by us- 
ing percolation theories Critical phenomena near the percolation thresh- 
old have been widely investigated in the electrical breakdown of granular metals or 
conductor-insulator composite The associated critical exponents have been 
,iQ OOEB0S0 and theoretically studied using continuum or lattice percolation 
deJ-li-i. In particular, large attention has been devoted to the critical expo- 
nents describing the resistance and its relative noise in terms of the medium properties 
(e.g. conducting particle fraction, defect concentration, etc. )0 — EI'El'i3 — El. However, very few 
attempt£i _ 0ii have been made so far to describe the behavior of a disordered medium over 
the full range of the applied stress, i.e. by studying the response of the system to an external 
bias when the bias strength covers the full range of linear and nonlinear regimes. Therefore, 
a satisfactory understanding of breakdown phenomena over the full dynamical regime is still 
missingEl'El. On the other hand, relevant information can be obtained from such a study, 
like: precursor phenomena, role of the disorder, existence of scaling laws, predictability of 
breakdown, etc.0-|Offl. 

Here we present a percolative model of sufficient generality to address the above issues. 
Our aim is to provide a theoretical framework to study response and fluctuation phenomena 
under linear and nonlinear regimes in a wide class of disordered systems. To this purpose, 
we study the evolution of a random resistor network (RRN) in which two competing pro- 
cesses are present, defect generation and defect recovery, which determine the values of the 
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elementary network resistances. Both processes are driven by the joint effect of an electri- 
cal bias and of the heat exchange between the network and the thermal bath. The bias is 
applied through a constant voltage or, alternatively, a constant current. Monte Carlo (MC) 
simulations are performed to investigate the network evolution in the full range of bias val- 
ues. Depending on the bias strength, an irreversible failure or a stationary state of the RRN 
can be achieved. By focusing on the steady-state, we analyze the behavior of the average 
network resistance, < R >, and the properties of the resistance fluctuations as a function of 
the bias. In constant voltage conditions, a scaling relation is found between < R > / < R >o 
and V/Vq, where < R >o is the linear regime resistance and Vq the threshold value for the 
onset of nonlinearity. A similar relation is found in constant current conditions. The relative 
variance of resistance fluctuations also exhibits a strong nonlinearity in the pre-breakdown 
regime whose properties are discussed. The fluctuation analysis is completed by a further 
investigation of the noise resistance spectrum and of the Gaussian features of the fluctuation 
amplitudes. Theoretical results agree with electrical breakdown measurements in thin films 
of compositesBli! and of other conduct ing0Jii@ or insulating materials!!. 

The paper is organized as follows. In Sect. 2 we briefly describe the model used. Section 
3 presents the results of the MC simulations for the resistance and its fluctuations. The 
main conclusions are drawn in Sect. 4. 



II. MODEL 

We study a two-dimensional random resistor network of total resistance R, made of N to t 
resistors, each of resistance r n , disposed on a square lattice. We take a square geometry, 
N x N, where N determines the linear size of the lattice, with the total number of resistors 
being N tot = 2N 2 . For the comparison with resistivity measurements in thin films, the value 
of N can be related to the ratio between the size of the sample and the average size of the 
grains composing the sample. An external bias, represented by a constant voltage V or by 
a constant current /, is applied to the RRN through electrical contacts realized by perfectly 
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conducting bars at the left and right hand sides of the network. A current i n is then flowing 
through each resistor. The RRN interacts with a thermal bath at temperature T and the 
resistances r n are taken to depend linearly on the local temperatures, T n , as: 

r n (T n ) =r [l + a(T n -T )] (1) 

In this expression a is the temperature coefficient of the resistance and T n is calculated by 
adopting a biased percolation modeH 



as: 



T n = T + A[r n e n + — - J2 M - r ^n)] (2) 

Here, N neig is the number of first neighbours around the nth resistor, the parameter A, 
measured in (K/W), describes the heat coupling of each resistor with the thermal bath 
and it determines the importance of Joule heating effects. The parameter B is taken to be 
equal to 3/4 to provide a uniform heating in the perfect network configuration. We note that 
Eq. (0) implies an instantaneous thermalization of each resistor at the value T n , therefore, by 
adopting Eq. (|^), for simplicity we are neglecting time dependent effects which are discussed 
in Ref. [ 0. 

In the initial state of the network (corresponding to the perfect network configuration 
with no heating) all the resistors are identical: r n = r^. Now, we assume that two competing 
percolative processes act to determine the RRN evolution. The first process consists of gen- 
erating fully insulating defects (resistors with very high resistance, i.e. broken resistors) with 
probabilityil'ii Wo = exp{—E D /K B T n ), where E D is an activation energy characteristic of 
the defect and Kb the Boltzmann constant. The second process consists of recovering the 
insulating defects with probability Wr = exp(— Z^fj/i^Tn), where Er is an activation en- 
ergy characteristic of this second process. For A^O, Eq. (0) implies that both the processes 
(defect generation and defect recovery) are correlated percolations. Indeed, the probability 
of breaking (recovering) a resistor is higher in the so called "hot spots" of the RRN0. On the 
other hand, for A = Eq. (Q) yields T n = T , which corresponds to random percolations@@. 
The same is true for vanishing small bias values, when Joule heating effects are negligible. 



As a result of the competition between these two percolations and depending on the 
parameters related to the particular physical system (E^, Er, A, a, r , N) and on the 
external conditions (specified by the bias conditions and the bath temperature), the RRN 
reaches a steady state or exhibits an irreversible breakdown. In the first case, the network 
resistance fluctuates around an average value < R >. In the second case, a critical fraction 
of defects p c , corresponding to the percolation threshold, is reached, i.e. R diverges due 
to the existence of at least one continuous path of defects between the upper and lower 
sides of the network^. We note that in the limit of a vanishing bias (random percolation) 
and infinite lattices (N — ► oo), the expression: Er < Ed + KbT ln[l + exp(— Ed/KbT )] 
provides a sufficient condition for the existence of a steady stated. 

The evolution of the RRN is obtained by MC simulations carried out according to the 
following iterative procedure, (i) Starting from the perfect lattice with given local currents, 
the local temperatures T n are calculated according to Eq. (H); (ii) the defects are generated 
with probability Wjj and the resistances of the unbroken resistors are changed as specified 
by Eq. ([]]); (iii) the currents i n are calculated by solving Kirchhoff's loop equations by the 
Gauss elimination method and the local temperatures are updated; (iv) the defects are 
recovered with probability Wr and the temperature dependence of unbroken resistors is 
again accounted for; (v) R, i n and T n are finally calculated and the procedure is iterated 
from (ii) until one of the two following possibilities is achieved. In the first, the percolation 
threshold is reached. In the second, the RRN attains a steady state; in this case the iteration 
runs long enough to allow a fluctuation analysis to be carried out. Each iteration step can 
be associated with an elementary time step on an appropriate time scale (to be calibrated 
with experiments). In this manner it is possible to represent the simulation of the resistance 
evolution over either a time or a frequency domain according to convenience. 

In the simulations, as reasonable values of the parameters, we have taken: N = 75, 
r = 1 (fi), a = 10~ 3 (K" 1 ), A = 5 x 10 5 (K/W), E D = 0.17 (eV), E R = 0.043 (eV), 
and T = 300 (K) if not stated otherwise. The values of the external bias range from 
0.05 < I < 2.8 (A) under constant current conditions, and from 0.05 < V < 3.5 (V) under 
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constant voltage conditions. 



III. RESULTS 

Figures 1 (a) and 1 (b) report a sampling of the resistance evolutions coming out from the 
simulations under steady state conditions for the cases of constant voltage (Fig. 1 (a)) and 
constant current (Fig. 1 (b)), respectively. At increasing the external bias, a remarkable 
increase of the average resistance and of the amplitude of fluctuations is evident in both 
cases while approaching the breakdown. This is reached by applying an external bias just 
above the highest values shown in these figures. Precisely, the RRN becomes unstable under 
constant voltage already by applying a voltage just above 3.0 (V) and, under constant 
current, by applying a current just above 2.1 (A). 

In the following, Figs. 2 to 4 will detail the behavior of the average resistance while Figs. 
5 to 10 will focus on the results of resistance fluctuations. 

A. Resistance 

Figure 2 reports the average value of the RRN resistance as a function of the applied 
bias (current or voltage). Each value is calculated by considering the time average on a 
single steady-state realization and then averaging over 20 independent realizations (RRNs 
subjected to the same bias conditions). The numerical uncertainty is found to be within 
0.01% at worst. At the lowest biases the resistance takes a value < R > which represents 
the intrinsic linear response property of the network (Ohmic regime)^. We note that for 
the activation energy values here chosen, the average fraction of defects at the lowest bias 
< p > is very small (< p > = 0.545 x 10 -2 ) so that < R >= 0.9978 (fi) is about 1% 
greater than the perfect network resistance. Under constant current conditions, when the 
current overcomes a certain value, Jo, the average resistance starts to become dependent on 
the bias. Thus, Iq sets the current scale value for the onset of nonlinearity. By adopting as 
criterium for the onset of nonlinearity a resistance increase of 0.05% over the linear response 
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value, we determine I = 0.090 ± 0.005 (A). For currents above Jo the resistance increase is 
smooth at moderate bias and it exhibits a sharp ramp (typical of a catastrophic behavior) at 
high bias, until a threshold current value lb = 2.10 (A) above which the RRN undergoes an 
irreversible breakdown. A step of current values 51 = 0.05 (A) has been used to determine 
lb, thus this value of SI gives the maximum uncertainty on lb- Accordingly, we have found 
h/Io = 23±1 and < R >b / < R >o= 1.6±0.1, where < R >b is the last stable value of the 
resistance calculated before the breakdown. We remark that the uncertainty on the ratio 
< R >b I < R >o is mainly due to the uncertainty on < R >b which reflects, amplified by 
the nonlinearity, the uncertainty on Similarly, under constant voltage conditions we have 
found V b /Vo = 32 ± 2 and < R > b / < R > = 1.5 ± 0.1, where V = 0.095 ± 0.005 (V) and 
Vb = 3.00 (V) are the voltage values corresponding respectively to the nonlinearity onset 
and to the electrical breakdown. The maximum uncertainty on Vb is 5V = 0.05 (V). The 
significantly higher value of the ratio Vb/Vo, when compared with that of the ratio h/h, 
is a quantitative indication that the system is more robust when biased under constant 
voltage than under constant current conditions. This property is further emphasized by the 
fact that the increase of the resistance in the pre-breakdown region exhibits a lower slope 
under constant voltage than under constant current conditions. It must be noticed that in 
spite of the significant difference of the ratios Vb/Vo and h/Io, the ratio < R >b / < R >o 
remains the same, within the error, under the different bias conditions. This result agrees 
with measurements of the ratio < R>b / < R >o performed in composites under the Joule 
regimeEl. 

To better analyze the dependence on the bias of the average resistance, Figs. 3 (a) and 3 
(b) report the log-log plot of the relative variation of the resistance as a function of V/Vq and 
I/Io, respectively. Figure 3 (a) shows that the relative variation of the average resistance 
scales with the ratio V/Vq in the whole region of applied voltages up to breakdown as: 



where a = (1.1 ± 0.1) x 10 4 is a dimensionless coefficient and 9 — 2.1 ± 0.1. Figure 3 (b) 




(3) 



S 



shows that also under constant current the relative variation of the average resistance scales 
with the ratio I/Iq in the moderate bias region, as: 

<R> - iwfir (4> 



< R > \h, 

with a' = (3.8±0.1) x 10 -4 and the same exponent 9 of Eq. (|3[), within numerical uncertainty. 
However, in this case we notice that in the pre-breakdown region the relative variation of the 
average resistance exhibits a superquadratic behavior characterized by a power law (I/Iq) 01 
with Qi = 3.7 ±0.1. 

This behavior of the average resistance with the applied bias can be understood as 
follows. From a macroscopic point of view, in a degradation process associated with a 
resistance increase, constant current conditions lead, at increasing bias, to a superquadratic 
increase of the power dissipated through Joule heating and thus to a major efficiency in the 
defect generation. By contrast, constant voltage conditions lead, at increasing bias, to a 
growth of the dissipated power which is subquadratic in the applied voltage, thus implying 
a minor efficiency in the generation of defects. In this respect, it must be underlined that 
the resistance increase driven by both kind of bias is further enhanced by the positive value 
of the temperature coefficient of the resistance. 

From a microscopic point of view, the RRN resistance depends on the average fraction 
of defects, < p >, according to the expression^: 

< R >~ | < p > -p c (5) 

where, for biased percolation, p c and p are functions of the bias strengthHH. Of course, for 
vanishing bias, p c and p take the well known values corresponding to random percolation: 
p c = 0.5 (for a square lattice) and p = 1.303 (universal value)0. On the other hand, the 
dependence of p c and p on the bias makes the analysis of < R > in terms of Eq. (H) quite 
problematic. Nevertheless, it is interesting to analyze the behavior of < p > on the bias 
strength and, to this purpose, Fig. 4 reports on a log-log plot the relative variation of < p > 
as a function of the normalized bias. As evidenced by the figure, for the case of constant 



voltage conditions the average fraction of defects scales with the same exponent 9 over the 
full range of bias values, as: 

<P>o \VoJ 

with b = (5.8 ± 0.3) x 10~ 3 a dimensionless coefficient. A similar expression holds under 
constant current conditions at moderate bias: 

^=l + tf(f)' (7) 
<P>o \IqJ 

with b' = (5.3 ±0.3) x 10~ 3 . However, under constant current, there is a significant increase 
of the slope in the pre-breakdown region up to a nearly above cubic power law. 

The behavior of the average fraction of defects in the moderate bias region can be un- 
derstood by generalising the results obtained for random stationary percolations!! to the 
presence of a moderate external bias responsible of Joule heating. In fact, for a RNN in 
a stationary state resulting from the competition of two random percolations, it has been 
foundS that < p >= x/(l + x), with x = W 0D (1 - W 0R )/W 0R , where W d and W 0R are 
the probabilities defined in section II with T n = Tq. At moderate bias, when the defect 
distribution is rather homogeneous and the variation of < R > small, we can take x = x(T) 
with T = Td + AT, where in the spirit of a mean-field theory, the average temperature 
increase can be expressed as AT = A < R > I 2 /(2N 2 ). Therefore, by differentiating: 



which by simple manipulation gives: 



A < p > 
<P>o 



A < R > H dx 



I 



2 



(9) 



_2iV 2 <p> [(l + x) 2 ] \dTj 

The term in the squared brackets of Eq. (|9|) can thus be identified with the V coefficient 
in Eq. ([?]). In fact, by introducing the parameters used in the simulations we find for the 
squared brackets in Eq. ([5|) the value 5.5 x 10~ 3 which well agrees with the previous reported 
value of b' . For constant voltage the factor A < R > Iq inside Eq. (|jD must be replaced 
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with AVq I < R >o- This gives for the square bracket in Eq. (^) the value of 6.1 x 10~ 3 , 
which again well agrees with the value of b in Eq. (||). Therefore, at moderate bias, the 
quadratic behavior of < p > is well explained by Eq. (|9[). In the pre-breakdown region, 
under constant current conditions, the significant increase of the resistance and the defect 
filamentation characteristic of biased percolationil'il add to give a superquadratic increase 
of < p > with the bias. By contrast, under constant voltage, these two effects act in an 
opposed manner, thus their compensation results in a quadratic dependence of < p > over 
the full range of bias values. 

From the above behavior of the average defect fraction we can now understand the bias 
dependence of < R >. At moderate bias (I/Io or V/V < 10) the increase of < R > is rather 
small. Accordingly, super and subquadratic effects on the dissipated power are negligible 
and the average fraction of defects over the intrinsic value grows quadratically with the bias 
strength (Eq. (|9])); at these bias values < p > <C 1 thus < R > oc < p > and also the av- 
erage resistance grows quadratically with bias. At high bias (i.e. I / 1 or V JVq > 10), when 
entering the pre-breakdown region, the increase of < R > is significant. Accordingly, super 
and subquadratic effects on the dissipated power are relevant as well as the increase of the 
average fraction of defects which now takes a nonuniform distribution (biased percolation). 
As a consequence, the percolation threshold is lowered with respect to the random perco- 
lation value!!. Depending on constant current or constant voltage conditions, this biased 
percolation effect adds to further enhancing the increase of the resistance well above the 
quadratic law or to keep the quadratic behavior, respectively. 

This interpretation is confirmed by MC simulations performed with the same parameters 
but setting the temperature coefficient of the resistance a = 0. In this case, we have found 
that the onset of nonlinearity as well as the electrical breakdown occur at higher values 
of the bias. Precisely we have found I' = 0.25 ± 0.03 (A), I' b = 2.65 ± 0.05 (A) (thus 
I'JI'q = 11 ± 1), < R > 6 / < R > = 1-2 ± 0.1 under constant current and V$ = 0.25 ± 0.03 
(V), V b ' = 3.25 ± 0.05 (V) (thus V b /V = 13 ± 1) and < R > b / < R > = 1.3 ± 0.1 under 
constant voltage conditions. The quadratic regime is always present at moderate bias values. 
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However, in the pre-breakdown regime we have found d\ = 3.1 ± 0.1 for constant current 
and 8y = 2.8 ± 0.1 for constant voltage conditions, respectively. 

We finally conclude this section by remarking that the behavior described by Eqs. (|3]) 
and (^) well agrees with measurements of the resistance of composites made in the Joule 
regime up to breakdown^. 



B. Fluctuations 

It is well known that important information about the properties and the stability of dif- 
ferent systems (physical, biological, social, etc.) can be obtained by studying the fluctuations 
around the average value of some characteristic quantities of the system0-!>0-§>§-§<E|-Ei 
Several features of the fluctuations are usually considered, and among these, the most im- 
portant are the variance (or the relative variance), the power spectrum and the Gaussianity 
property. Here, we analyze the fluctuations of the RNN resistance in the full range of bias 
values, devoting particular attention to the features that can be identified as failure precur- 
sors. Figure 5 reports the relative variance of resistance fluctuations, £ =< 5R 2 > / < R > 2 , 
as a function of the applied bias. The data correspond to the same simulations presented 
in Fig. 2. The same procedure of time averaging over a single simulation and then making 
ensemble averages over 20 realizations provides an uncertainty of 3% at worst. At low bias, 
the relative variance is found to achieve a constant value S = (4.0 ± 0.1) x 10~ 6 . This 
value represents an intrinsic property of the system determined by the competition of two 
random percolations, as already discussed for the case of the resistance. A detailed anal- 



ysis of the scaling properties of S can be found in Ref. [ |3~9 |. At increasing bias, when 
/ > Iq or V > Vq, the systematic increase of £ reveals the existence of a nonlinear regime 
of the response. Similarly to the behavior of the resistance, by approaching the electrical 
breakdown E exhibits a significant increase which is steeper under constant current than 
under constant voltage conditions. This feature is ascribed to the better stability of the 
RRN under constant voltage, since resistance fluctuations in excess over the mean value are 
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damped in this condition. By contrast, the same kind of fluctuations are enhanced under 
constant current conditions. To detail the dependence of £ on the bias, we have reported in 
Fig. 6 (a) its relative variation, [£y — E ]/E , as a function of V/Vq under constant voltage 
conditions. Two regions can be identified in the nonlinear regime: a moderate bias region 
(V/Vo < 10), characterized by a quadratic dependence on the applied voltage, and a pre- 
breakdown region at the highest voltages (V/Vo > 10), where a superquadratic dependence 
up to a cubic power law is evidenced. Accordingly, in the first region it is: 

with c = (1.14 ±0.06) x 10~ 2 a dimensionless coefficient, while in the second region we have 
found: 

Y!„ /V\ f lv 

(11) 



S ~ Wo, 

where rj v = 3.0 ± 0.1. This behavior of E under constant voltage conditions should be 
compared with that of the average resistance which remains always a quadratic function up 
to the breakdown, as shown in Fig. 3 (a). The emergence of a superquadratic behavior of 
E in the pre-breakdown region reflects the higher sensitivity of the resistance fluctuations 
with respect to the average value of the resistance to the instability of the network. The 
same analysis performed for the case of constant current is shown in Fig. 6 (b). Again, at 
moderate bias [Ej — E ]/E increases quadratically with the current as: 



Eo \Io 

with d = (6.6±0.4) x 10~ 3 a dimensionless coefficient. On the other hand, the superquadratic 
dependence characterizes the pre-breakdown region as: 



E \/o/ 

with r)i = 5.4 ± 0.1. We note, that rji is significantly greater than rjv, according to the 
behaviors shown in Fig. 5. Moreover, r\i is greater than the resistance exponent in the 
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pre-breakdown region dj. It is noteworthy, that by neglecting the effect of the temperature 
coefficient of the resistance, i.e. by taking a = 0, the moderate bias region remains charac- 
terized by a a quadratic increase which is common to both constant current and constant 
voltage conditions. However, in the pre-breakdown region we have found r)v = rjj = 4.0 ±0.1. 
Thus, the above results prove that the quadratic dependence of £ on the bias in the mod- 
erate bias region is a feature independent of the conditions on which the bias is applied and 
also independent of the value of the thermal coefficient of the resistance. This reflects the 
fact that at moderate bias the two following conditions are satisfied: i) the variation of the 
average fraction of defects follows Eq. @; ii) the average fraction of defects is much smaller 
than the percolation threshold value, i.e. < p > <C p c . Therefore, in this nearly perfect 
regime S oc < p > oc < R >B On the other hand, the superquadratic behavior of £ in the 
pre-breakdown region is a consequence of the fact that at high biases the local correlations 
in the defect generation and recovery processes become strong. Defect filament at ion, char- 
acteristic of biased percolationil'il, implies that a small variation in the number of defects 
bring strong resistance fluctuations. 

Figure 7 reports the relative variance of resistance fluctuations normalized to the linear 
regime value as a function of the normalized resistance. We can see that the former quantity 
exhibits an increase over nearly two orders of magnitude in the same range of bias where the 
latter quantity increases for about 60% only. The log-log plot of [X — S ]/S as a function of 
[< R > — < R > ]/ < R >o, shown in Fig. 8, confirms the proportionality of S with < R > 
in the moderate bias region and evidences the following power-law in the pre-breakdown 
region. 

g ~ So] ^ [<R>-<R> ] <: 
S < R >o 

where the value of the exponent, ( = 1.5 ±0.1, within the statistical uncertainty is the same 
for constant voltage or constant current conditions. This result is consistent with Eqs. (|3]) 
and (|). 

To complete the study of resistance fluctuations, we have investigated the Gaussian 
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properties of the fluctuation amplitudes and the spectra in the frequency domain. Figures 9 
(a) and 9 (b) report the distribution function of the resistance fluctuations p(R) for different 
currents and different bath temperatures, respectively. The dotted and the long-dashed 
curves in Fig. 9 (a) correspond to fit with a Gaussian distribution the data corresponding to 
different current values in the pre-breakdown region at T = 300 (K), while in Fig. 9 (b) to 
fit the data corresponding to different bath temperatures at / = 1.5 (A). When approaching 
the breakdown conditions, at increasing current in Fig. 9 (a) and at increasing temperature 
in Fig. 9 (b), the simulations show the onset of a non-Gaussian behavior characterized by 
the enhancement of the probability for the positive fluctuations with respect to the Gaussian 
distribution. The emergence of a non-Gaussian behavior near the breakdown, in agreement 
with experiments0'0'il, can be considered a relevant precursor of failure. 

Figures 10 (a) and 10 (b) show the spectral densities of resistance fluctuations for the 
same conditions of Fig. 9 (a) and 9 (b). The spectral densities have been calculated 
by Fourier transforming the corresponding correlation functions C$n(t). We have found 
Lorentzian spectra in all the cases. For a given temperature, within the numerical un- 
certainty, the corner frequency is found to be independent of the applied current, while the 
value of the plateau increases more than quadratically with the current in the pre-breakdown 
region. For a given current in the pre-breakdown region, both the corner frequency and the 
value of the plateau are found to increase at increasing temperatures. This fact indicates 
that the characteristic times of fluctuations while depending on temperature they are inde- 
pendent of the applied bias. We finally remark that the power spectral density in Figs. 10 
(a) and 10 (b) are in good agreement with experiment 



IV. CONCLUSIONS 

We have studied the stationary state of a two-dimensional RRN resulting from the compe- 
tition between two biased percolations. The two percolative processes consist of the breaking 
and recovering of elementary resistors and they are driven by the joint effects of an electrical 
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bias and of the heat exchange with a thermal bath. The electrical bias is set up by applying 
a constant voltage or a constant current. MC simulations have been performed to analyze 
the network resistance and its fluctuation properties in the full range of bias values, covering 
linear and nonlinear regimes up to the breakdown limit. The nonlinear regime starts for 
biases greater than the threshold values, Vq or J (nonlinearity onset values), and it extends 
until the values Vb or lb (breakdown values). The ratios Vb/Vo, h/Io, < R >b / < R >o can 
be considered as relevant indicators characterizing the breakdown properties of the system. 
We have found that the ratio < R >b / < R >o is independent of the bias conditions but it 
depends on the thermal coefficient of the resistance. This result agrees with measurements of 
this ratio performed in compositesEl. Moreover, we have found that under constant voltage 
conditions the relative variation of the average resistance scales quadratically with the ratio 
V/Vq over the full nonlinear regime. A similar relation has been found under constant cur- 
rent conditions, but in this case a superquadratic dependence emerges in the pre-breakdown 
region (I > 10 Jo)- F° r what concerns the relative variance of resistance fluctuations we 
have found that at moderate bias (Iq < I < 10 Iq) it grows quadratically with the bias, 
independently of constant-current or constant-voltage conditions. On the other hand a su- 
perquadratic dependence appears in the pre-breakdown region. The presence of two distinct 
regions in the nonlinear regime has been explained in the following terms. At moderate bias 
there is a small Joule heating giving rise to the competitions of two nearly random percola- 
tions and to a small variation of the RRN resistance. At high bias, local correlations typical 
of biased percolation become important. Moreover, also the resistance variation becomes 
strong. Therefore, as a consequence of both these effects the resistance noise depends su- 
perquadratically on the bias. The power spectral density of resistance fluctuations is found 
to exhibit a Lorentzian spectrum and in the pre-breakdown region the amplitude of fluctua- 
tions shows the onset of non-Gaussian behaviors. Theoretical results agree qualitatively, and 
in some cases quantitatively, with breakdown experiments performed in compositesBS and 
in other conductingSBIl or nonconducting materialsil. We conclude that the study of the 
stationary state of a RRN, resulting from the competition of biased percolations, provides 
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an unifying framework for the interpretation of several non-linear transport phenomena in 
a variety of disordered systems. 
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FIGURE CAPTIONS 



Fig. 1 (a) - Resistance evolutions of a steady state 75 x 75 network under constant voltage 
conditions for increasing values of the bias. Going from the lower to the upper curves, the 
voltage is V = 0.2, 0.5, 1.0, 1.5, 2.0, 3.0 (V). The values of all the other parameters are 
specified in the text. 

Fig. 1 (b) - Resistance evolutions of a steady state 75 x 75 network under constant current 
conditions for increasing values of the bias. Going from the lower to the upper curves, the 
current is: / = 0.2, 0.5, 1.0, 1.5, 2.1 (A). The values of all the other parameters are specified 
in the text. 

Fig. 2 - Average network resistance as a function of the external bias. The averages are 
calculated over 20 networks of sizes 75 x 75 subjected to the same bias conditions. The circles 
correspond to biasing the network by an external constant current, the triangles correspond 
to a constant voltage bias. 

Fig. 3 (a) - Log-Log plot of the relative variation of the network resistance calculated under 
constant voltage as a function of the ratio V/V . The data are the same of those in Fig. 2 
and only the nonlinear regime is shown. The dotted line represents the power-law fit with 
a value for the exponent 9 = 2.1 ±0.1. 

Fig. 3 (b) - Log-Log plot of the relative variation of the network resistance calculated under 
constant current as a function of the ratio / / 1 . The data are the same of those in Fig. 2 
and only the nonlinear regime is shown. The long-dashed line represents the power-law fit 
of the data in the moderate bias region, the value of the exponent is 9 = 2.1 ±0.1. The solid 
line fits the data in the pre-breakdown region with a power-law, the value of the exponent 
is 6i = 3.7 ±0.1. 

Fig. 4 - Log-Log plot of the relative variation of the average number of defects under constant 
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voltage (triangles) and under constant current conditions (open circles) as a function of the 
normalized bias (V/V or I/Iq). The solid line represents the power-law fit with a value 
for the exponent 6 = 2.1 ± 0.1. The long-dashed line fits constant-current data in the 
pre-breakdown region with a power-law of exponent 3.3 ±0.1. 

Fig. 5 - Relative variance of resistance fluctuations, £, under constant voltage conditions 
(triangles) and under constant current conditions (open circles) as a function of the external 
bias. Each point has been obtained by averaging the variances calculated for 20 networks 
of sizes 75 x 75 subjected to the same bias conditions. 

Fig. 6 (a) - Log-Log plot of the relative variation of £ under constant voltage conditions 
as a function of V/V . The long-dashed line represents the power-law fit to the data in the 
moderate bias region, the value of the exponent is 2.0 ±0.1. The solid line fits data in the 
pre-breakdown region with a power-law of exponent 3.0 ± 0.1. 

Fig. 6 (b) - Log-Log plot of the relative variation of £ under constant current conditions 
as a function of I / Iq. The long-dashed line represents the power-law fit to the data in the 
moderate bias region, the value of the exponent is 2.1 ±0.1. The solid line fits data in the 
pre-breakdown region with a power-law of exponent 5.4 ±0.1. 

Fig. 7 - Relative variance of resistance fluctuations normalized to the linear regime value 
as a function of the normalized average resistance. Open circles refer to constant current 
condition and triangles to constant voltage condition, respectively. 

Fig. 8 - Log-Log plot of the relative variation of X as a function of the relative variation of 
the average resistance under constant voltage condition (open circles) and constant voltage 
condition (triangles). The long-dashed line gives the best fit with a power-law with exponent 
1.0 ± 0.1. The solid and the dotted lines fit with a power-law the data in the pre-breakdown 
region, the exponent are 1.5 ± 0.1 and 1.6 ± 0.1, respectively. 
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Fig. 9 (a) - Distribution function of the resistance fluctuations for two values of the applied 
current. Open circles refer to I = 1.5 (A) and full diamonds to / = 1.8 (A). The scale is 
a linear-log, therefore the dashed and the long-dashed curves, coming from parabolic fits of 
the two sets of data, correspond to Gaussian distributions. 

Fig. 9 (b) - Distribution function of the resistance fluctuations under constant current 
conditions for two values of the substrate temperature. Open circles refer to T = 300 (K) 
(the same data reported in Fig. 9 (a) ) and full triangles to T = 450 (K). The scale is 
linear-log, therefore the dashed and the long-dashed curves, coming from parabolic fits of 
the two sets of data, correspond to Gaussian distributions. 

Fig. 10 (a) - Power spectral density of resistance fluctuations under constant current con- 
ditions for / = 1.5 (A) (dotted curve) and for / = 1.8 (A) (solid curve). In both cases the 
substrate temperature is 300 (K). 

Fig. 10 (b) - Power spectral density of resistance fluctuations under constant current con- 
dition for T = 300 K (dotted curve) and for T = 450 (K) (solid curve). In both cases the 
bias current is 1.5 (A). 
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